Competitive nucleation and the Ostwald rule 
in a generalized Potts model with multiple metastable phases 
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We introduce a simple nearest-neighbor spin model with multiple metastable phases, the number and decay 
pathways of which are explicitly controlled by the parameters of the system. With this model we can construct, 
for example, a system which evolves through an arbitrarily long succession of metastable phases. We also 
construct systems in which different phases may nucleate competitively from a single initial phase. For such a 
system, we present a general method to extract from numerical simulations the individual nucleation rates of the 
nucleating phases. The results show that the Ostwald rule, which predicts which phase will nucleate, must be 
modified probabilistically when the new phases are almost equally stable. Finally, we show that the nucleation 
rate of a phase depends, among other things, on the number of other phases accessible from it. 

PACS numbers: 64.60.My, 64.60.Qb, 05.10.Ln, 05.50.+q 



Metastability is a ubiquitous phenomenon in nature. 
Broadly speaking, it occurs when a system is "trapped" in 
a phase different from equilibrium. This non-equilibrium 
phase, the metastable state, can last for extremely long times. 
Thus, it is not surprising that metastable states play a cru- 
cial role in many physical processes and are at the center of 
much current research. For example, recently an intermedi- 
ate metastable phase was shown to provide an easier pathway 
for the growth of crystal nuclei from fluids {nucleation), with 
implications for the crystallization of proteins iQllil]. Proteins 
themselves are known to get stuck in misfolded metastable 
structures |3j|, preventing them from reaching their equilib- 
rium configuration. The phenomenology observed in these 
and in many other systems can be thought of as arising from a 
complicated energy landscape, with several local "metastable 
minima" where the trapping occurs J3l- The extreme situation 
is that of glasses, in which the energy landscape can have ex- 
temely many local minima hindering relaxation of the system 
to a thermodynamically stable crystal ^ . 

The above systems present at least several metastable states. 
These states and the transitions between them usually arise 
from the microscopic interactions in a complicated way. 
When this is the case, the study of phenomena such as compe- 
tition between nucleating phases and specific nucleation path- 
ways may be obscured. In view of this, in this work we 
present a simple spin model with nearest-neighbour interac- 
tions, where the number of metastable phases and the decay 
pathways between them can be explicitly specified by varying 
the model parameters. It thus serves as a test-bed for theoreti- 
cal results relating to systems with multiple metastable phases 
ifiHIilL just as the kinetic Ising model, a special case of our 
model, has been central in the study of systems with a sin- 
gle metastable phase [9]. As discussed below, the model also 
describes the adsorption of multiple chemical species onto a 
surface, an interesting physical problem in its own right. 

After presenting the model, as an illustration of a possible 
application, we construct a system with arbitrarily long suc- 
cessions of metastable states. We then focus on competition 
between phases nucleating from a single initial metastable 



phase. An important question in this context is to understand 
which phases nucleate under which conditions. The Ostwald 
rule states that the nucleating phase is the one with the small- 
est free energy barrier from the initial phase: see Ref. [10] 
and references therein. Previous results have supported this 
prediction flUl . 

We show that in general the Ostwald rule must be modified 
probabilistically when the new phases are of similar stabil- 
ity, using an argument based on individual nucleation rates of 
each phase. We give a method by which these rates can be 
measured in simulations or experiments, and show that there 
is a parameter regime in which any of the new phases may 
nucleate — only the nucleation probability of each phase can 
be established, with the outcome in any given run being un- 
predictable. We finally show that the nucleation probability of 
a phase depends on the phases accessible from it. 

Model details:- Our model is based on the Potts model, in 
which each spin has one of q states [12] and each phase has a 
majority of spins in one state; the Ising model corresponds to 
q = 2. The relative stability of each phase is controlled by ex- 
ternal fields, and the interplay of these fields with interactions 
between different spin states allows us to obtain any desired 
transition pathways between phases. 

Viewing the fields as chemical potentials, we can recast the 
model as a multi-component lattice gas which describes ad- 
sorption on a lattice substrate (e.g., a crystal plane) of mul- 
tiple chemical species with lateral interactions II 1311 . Much 
experimental work has been done on the thermodynamics of 
such systems, but little on the kinetics — see |14l| and refer- 
ences therein; nonetheless, our results should be testable in 
that context. A more complicated system where the kinetics 
has been characterized is a colloid-polymer system 1 15, Trill , 
where possible pathways were found from considerations of 
the free energy landscape [17]. Our approach is complemen- 
tary in that specific pathways result from microscopic interac- 
tions. 

We work on an L x L square lattice with N := L 2 spins and 
periodic boundary conditions, although the results are qualiti- 
tavely unaffected by lattice type. Each lattice site i has a spin 
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Oj taking values in {1, . . . ,q}, and the energy of a configura- 
tion 6 is given by the Hamiltonian 



H(&) 



ce=l 



(1) 



Here, M a := £ ; a is the magnetisation (^number of spins) 
of the spin type a; 8 a ,y = 1 if 0£ = 7, and otherwise. The 



first term is a sum over nearest-neighbor pairs of spins of a 
symmetric interaction energy J a y — Jy. a , and the second de- 
scribes the effect of external fields h a acting on spin type a. 

We set the diagonal elements J a a of the interaction ma- 
trix to unity (Ja.a '■= 1 for all a), so that in the absence of 
non-diagonal interactions and fields, the model reduces to the 
standard Potts model 111 211 . This has q symmetrical phases 
coexisting below a critical temperature T c = 1/ (ln(l + ^fq)); 
each phase has a majority of spins in one of the q spin states. 
Including fields breaks the symmetry between phases. If 
h a = hp then the a and 7 phases coexist, with a first-order 
phase transition between them, for T < T c ; this is at the ori- 
gin of metastability. Weak non-diagonal interactions do not 
qualititavely affect this coexistence. 

To evolve the system we choose discrete-time Metropolis 
dynamics 111 811 : at each time step, a spin and its new value 
are chosen at random, the increment AH of the Hamiltonian 
([T} for this change is calculated, and the update is accepted 
with probability min{l,exp(— fiAH)}, where J3 := 1/T is the 
inverse temperature. This gives a Markov chain on the space 
of all possible configurations. 

This Markov chain has a unique equilibrium distribution, 
concentrated on the phase(s) with the largest h a . The other 
phases are metastable, that is, when started in such a phase a, 
the system stays there for some time, before a transition to a 
more stable phase 7 is nucleated by the appearance of a crit- 
ical droplet of the 7 phase. At sufficiently low temperatures, 
the relative stability is determined by hy > h a . The reverse 
transition is exponentially unlikely. 

In the standard Potts model, the equilibrium phase (almost) 
always nucleates. To obtain non-trivial transition pathways, 
nucleation of other phases must be promoted. This we achieve 
using non-diagonal interactions between distinct spin types 
a 7^ 7: setting J a .y > favors nucleation of 7 droplets inside 
the a phase by lowering the surface tension between a and 7 
regions, and hence decreasing the droplet free energy of for- 
mation {nucleation barrier), whereas formation of 7 droplets 
in the a phase is suppressed if J a y < 0. 

We can now construct models whose phases obey arbitrary 
metastable transition graphs. These are directed graphs with 
the restriction that no loops returning to a previously visited 
phase are allowed. Each vertex corresponds to one phase, la- 
beled by its dominant spin state, and each arrow to a desired 
transition: a — > 7 means that phase 7 can nucleate directly 
from phase a. Fig.Q]shows example transition graphs. 

To construct a model corresponding to a given transition 
graph, we proceed as follows. The number of spin types, q, is 
the number of vertices in the graph. To each spin type a we 
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FIG. 1: Metastable transition graphs: (a) kinetic Ising model; (b) suc- 
cession of 3 phases; (c) single metastable phase decaying to two com- 
peting phases; (d) as in (c), but such that one phase can decay further; 
(e) three competing phases. 
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FIG. 2: (Color online) Time dependence of magnetisation per site 
M a /N of each phase a, in a single run of the model with transition 
graph 1^2^3^4^5(a succession of phases). Parameters are 
L = 50, j8 = 1.25, h a = 0.1(a — 1), K\ =0.1, and K 2 = 1.0. Labels 
denote the dominant phase. Configuration snapshots depict post- 
critical nuclei of each new phase embedded in the previous phase. 
These grow to fill the system, producing the next phase in sequence. 



assign a field h a , with hy > h a if 7 is below a in the graph. 
The off-diagonal interactions are given by J a y : = K\ > (at- 
tractive) if a — > 7, and J a y := —K2 < (repulsive) otherwise. 
K2 must be large enough to inhibit immediate formation of 
non-adjacent phases with large fields. 

As an illustration, we construct a model exhibiting a lin- 
ear succession of metastable phases with transition graph 
1 — > 2 — > • • • — > q. We impose fields = h\ < h% < ■■■ < h q 
and attractive interactions J a ,a±\ '■= Ki > between neigh- 
bouring states, and set all other non-diagonal interactions to 
— K2 < 0. With suitable, moderately robust, parameters, we 
observe the desired behavior, shown for q = 5 in Fig. [2] 

A three-phase succession was previously observed in a ki- 
netic Blume-Capel model llSl 12011 . corresponding to a special 
case of our model with q = 3 112 ill . The physical reason for the 
observed transitions is, however, much more transparent with 
the Hamiltonian in the form ([T), with its intuitive interpreta- 
tion in terms of attractive and repulsive interactions. 

Competitive nucleation:- We now turn to the decay of one 
metastable phase into two competing phases (Fig.QJc)). Sear 
studied competitive heterogeneous nucleation (occurring on 
impurities) in the 3-state standard Potts model 111 111 . In con- 
trast, all behaviors discussed in this work are endogenous: ob- 
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FIG. 3: (Color online) Probability piQij) that phase 2 nucleates be- 
fore phase 3 as a function of Ah/li2, with hi = 0.1, K\ = 0.1 and 
K2 — 1, and jS = 1.25, L = 50 unless otherwise noted. Up to n = 10 4 
trials were used for each data point; statistical errors are of the or- 
der of the symbol size. Inset: system-size dependence of P2 for two 
values of with fixed parameter values. 
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FIG. 4: (Color online) (a) Comparison of nucleation rates A/ (A3) of 
phases i = 2, 3 from phase 1 in the model of Fig.fTJc) with I12 fixed, 
calculated directly from simulations and using forward-flux sam- 
pling, (b) Nucleation probability p 2 {h 2 ) for /13 = 0.1 and I14 = 0.2 
varying I12, with and without phase 4, in the model of Fig. [TJd). 
Dashed lines show positions of equal nucleation probability and 
equal field of the two phases. In both subfigures parameters are as in 
Fig.fJ with L = 50 and j8 = 1 .25. 



served transitions are not caused by external influences, but 
rather arise spontaneously from within the system itself. 

We fix = h\ < h 2 < /13, /i,2 = /i,3 > and J 2) 3 < 0. Let 
Ah := /i3 — h 2 be the field difference between the new phases, 
2 and 3. When Ah = 0, these phases are symmetrical, each 
nucleating half of the time, while for Ah > 0, we expect the 
1-3 free energy barrier to be lower than the 1-2 one, so that 
according to the Ostwald rule, only phase 3 should nucleate. 
To test this, we perform n simulations starting from phase 1 
for each Ah, in n 2 of which phase 2 nucleates before phase 3. 
The ratio n 2 /n is then an estimate of the probability p 2 (Ah) 
that phase 2 nucleates first. For efficiency' we use a rejection- 
free version of the Metropolis method R2U 12211 . 

Fig. [3]plots p 2 as a function of a non-dimensionalised Ah. 
For Ah sufficiently close to 0, phase 2 can still nucleate first, 
contrary to the simple Ostwald rule. The probability that it 
does so rapidly decreases for larger Ah, until a point beyond 
which phase 3 effectively always nucleates. 

To explain these results in a general context, we assume, as 
in classical nucleation theory M23I1 . that there are well-defined 
nucleation rates A;(L) of phases i = 2, 3, giving the number of 
critical nuclei which form per unit time in a system of size L. 
The nucleation rates per site are jU,-(L) := A,-(L) /N. 

A nucleation rate is the inverse of a mean nucleation time, 
which can be measured in experiments or simulations by aver- 
aging over many nucleation events in independent runs. In the 
case of competitive nucleation, however, we can only measure 
the mean time T for the first phase to nucleate, after which 
this phase invades the entire system. The rate of this first 
nucleation is X 2 + A3, since the total number of nucleation 
events per unit time is the sum of those for each type, so that 
T = 1/ (A-2 + A3). For convenience, in simulations 1 is taken 
to be the time until the new phase occupies half the system. 

Under the same assumptions, the probability that phase 2 



nucleates first is p 2 — Prob{72 < T3}, where 7), the time for 
phase i to nucleate, is an exponentially distributed random 
variable with mean 1/A,\ This gives p 2 = X 2 j (A2 + A3) = A2T. 

Individual nucleation rates of the two phases, which can- 
not be obtained directly, can now be calculated as: A2 = 
p 2 jx and A3 = (1 — p 2 )/T. This generalizes to P compet- 
ing phases, where the measurable quantities are the mean nu- 
cleation time T = 1/L^li A,- and the nucleation probabilities 
Pi = Aj/£/A/ = Mi/L/M; of each phase i. The nucleation 
rate of phase i is then A,- = pi/ x. 

Figure|4f a) shows A2 and A3 calculated in this way. To con- 
firm the validity of such calculations, we use the "forward-flux 
sampling" method B24 12511 . which directly calculates the tran- 
sition rate between two phases in a stochastic system. This 
has previously been used to study nucleation rates in the Ising 
model lErjl I27I1 . In our case, the possibility of escape to an 
additional new phase must be taken into account [21]. Fig- 
ure St a) shows that the results indeed coincide with those of 
the direct method, within statistical errors. We remark that we 
are unaware of any analytical prediction giving the observed 
variation of nucleation rates. 

The above considerations used "in reverse" confirm that 
the Ostwald rule must in general be modified when the new 
phases have similar stabilities, as follows. Consider phases 
which are equally stable for given parameter values. We ex- 
pect nucleation barriers, and hence nucleation rates, to vary 
continuously with the parameters, so that the nucleation prob- 
abilities pi also vary continuously. Hence there is a region, 
where the phases have similar stabilities, in which all nucle- 
ation probabilities are non-zero — only the probability of each 
phase nucleating is well-defined, with the outcome in any 
given run being stochastic, as in Fig. [3] The definite prediction 
given by the Ostwald rule is thus invalid in this region. 

To see how our results depend on system size L, we note 
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FIG. 5: (Color online) Configuration snapshots of the multidroplet 
regime in the model with transition graph Fig. [TJe), with L = 200, 
j8 = 1.25, and h\ = 0, ho, = I13 = h^ = 0.15. Nucleation of droplets 
of three equally stable phases from a single metastable phase is fol- 
lowed by droplet growth and then domain coarsening. The "coating" 
of the initial phase visible between domains in the final snapshot is 
due to repulsive interactions between the new phases. 

that in a broad region of L, the per-site nucleation rates 
and hence also the pi, are independent of L, as confirmed by 
the plateaus in the inset of Fig. [3] This is valid when the nu- 
cleation process is mediated by growth of a single droplet [9]. 
Note that this regime may be of relevance for macroscopically 
large systems |9J]. 

Above a certain system size, however, droplets of differ- 
ent phases may nucleate before any dominates the system 
(the 'multidroplet' regime) |U. This results in coarsening, 
as shown in Fig. [5] for three competing phases of equal sta- 
bility. Even if a phase-a droplet nucleates first, droplets of 
a more-stable phase may then nucleate and grow to dominate 
the system before the a phase can do so, thus reducing p a , as 
seen in Fig. [3] 

Finally, a generic non-symmetric case can be obtained by 
adding a new phase, 4, and a decay path, 3 — > 4, as in Fig.Qld). 
Even when /12 = ^3, phase 3 now nucleates more often, as 
shown by the horizontally displaced nucleation probability 
curve in Fig. Hfb): the presence of phase 4 reduces P2 by 
roughly half. This is due to an entropic effect: there are more 
3-dominated critical droplets than 2-dominated ones, since 
spins of type 4 also appear in the former droplets, resulting 
in nucleation of a binary mixture ll23ll : a similar effect is vis- 
ible in Fig. [TJ There is thus a lower free energy barrier to 
form phase 3, and yet nucleation of both phases 2 and 3 is still 
observed, again at odds with the simple Ostwald rule. 

In summary, we have introduced a generalized Potts model 
which can easily be tuned to have any given number of 
metastable phases and arbitrary transitions between them. We 
have shown generally that individual nucleation rates of com- 
petitively nucleating phases can be calculated from exper- 
imentally measurable quantities, and that the Ostwald rule 
must be modified when the nucleating phases have compa- 



rable stabilities. In future work [21], we will study the model 
in detail and compare its properties with theoretical results [8] 
on systems with multiple metastable phases. 
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